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CsJ ; ABSTRACT 

. Context. The true mass-loss rates from massive stars are important for many branches of astrophysics. For the correct modeling of the 

rj • resonance lines, which are among the key diagnostics of stellar mass-loss, the stellar wind clumping turned out to be very important. 

'"^ ■ In order to incorporate clumping into radiative transfer calculation, 3-D models are required. Various properties of the clumps may 

T-H , have strong impact on the resonance line formation and, therefore, on the determination of empirical mass-loss rates. 

^Nj . Aims. We incorporate the 3-D nature of the stellar wind clumping into radiative transfer calculations and investigate how different 

model parameters influence the resonance line formation. 

I I' Methods. We develop a full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. The number density 

I of clumps follows the mass conservation. For the first time, realistic 3-D models that describe the dense as well as the tenuous wind 

■ components are used to model the formation of resonance lines in a clumped stellar wind. At the same time, non-monotonic velocity 
• ■ fields are accounted for. 

■ Results. The 3-D density and velocity wind inhomogeneities show very strong impact on the resonance line formation. The different 
parameters describing the clumping and the velocity field results in different line strengths and profiles. We present a set of represen- 

Q ■ tative models for various sets of model parameters and investigate how the resonance lines are affected. Our 3-D models show that 

5_j ■ the line opacity is reduced for larger clump separation and for more shallow velocity gradients within the clumps. 

■ Conclusions. Our new model demonstrates that to obtain empirically correct mass-loss rates from the UV resonance lines, the wind 
I clumping and its 3-D nature must be taken into account. 

I Key words, stars: winds, outflows, clumping - stars: mass-loss - stars: early-type 



. 1 . Introduction cal simulations of line-driven stellar winds that have been per- 

■ formed using the simpl ifying assumption of spherical symme- 
; Hot massive stars lose mass via stellar winds. Mass loss plays j.j^ geometry (Feldmeier;T99?;'Feldmeieretall|lMi 
. an important role in the stellar evolutionand strongly affect the Q^Q^ki et al., 1988; Runacres & Owocki, 200^ or pseudo-2D 

^ . interstellar environment (| Bresolin etaU |2008|). In order to better geometry (iDessart & OwockiL l2003Ll2005i) . 

■ understand the evolution of hot massive stars and their influence 

(N| 1 on the environment it is necessai-y to determine their mass-loss Theoretical predictions ai-e supported by direct observational 

^ ■ rates reliably. Many sophisti cated models of the stella r evolu- evidence of clumping. Eversberg et al. ( 1998) found stochastic 

ILJ ■ tion have been calculated (see iMaeder & Mevnetl iSl for a re- variable structures in the He emission hne of ^ Puppis, which 

■ view), however relying on mass-loss rates, that are still far from drift with time from the line centre to the blue edge of the line, 
being empirically established These structures were explained by accelerated wind inhomo- 

Since the 1970s, significant progress in understanding geneities moving outwards. High-resolution spectroscopic mon- 

5^ ■ the physics of hot sta r winds has been made (see e.g . itoring of the line-profile variations (LPVs) in the emission lines 

Lucv & Solomoa iMdc ICastor et al.L [iffl iPauldrachet all °f "i"^ Wolf-Rayet (WR) stars dLepine & Moffatt [1999|), inves- 
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1986'), assuming the standard wind model (stationai-y, spheri- ligation of the LPVs of Ha for a large sample of O-type super- 

cally symmetric wind with uniform flow). Winds of hot stars gi^^^ (|Markova et al.U 2005), and direct spectroscopic observa- 

are driven by radiation absorbed in spectral lines, so-called the Hon of five O-type massive star of different evolutionary stages 

line-driven winds. For recent reviews of the physical mecha- <|Lepme & Mofi^att |200S), suggest that clumpy structures are a 

nism of hne driv ing an d of hne-dr iven stell ar winds see, e.g., common property and a universal phenomenon of all hot stai- 
iKrticka&KuMl (l2007h . |Pulsetan (l2U58h . or lOwockil (l2010i) . 

In the last decades there was a growing e videnc e that the In addition to the stochastic small-scale wind structures, 

stellar winds are not smooth dHamann et al. I l2008l) . Detailed there is also strong evidence for the presence of large-scale 

theoretical studies showed that the line-driven winds are in- wind structures, namely the discrete absorption components 

trinsically unstable (Lucv & White, 1980). Due to this instabil- (DACs). These DACs are observed to propagate bluewards 

ity shocks and wind density structures (clumps) develop in the through the UV resonance line profiles of nearly all O-type stars 

winds. Theoretical evidence of clumping is based on numeri- (iPrinia & HowarthL Il986t iHamann et all 1200 ih . To explain the 
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observed DAC properties qualitatively, the mod el of Co r otatin g 
Interaction Regions (CIR) has been proposed (iMuUanL Il984l) . 
The CIRs form in a rotating stars when high-density, low-speed 
win d streams colhde with low- density, high-speed streams (see, 
e.g.. lCranmer & Owockil[T996l) . 

One of the indirect evidences of clumping in hot star 
winds comes from X-rays observations. Detailed investiga- 
tion of the X-ray transport in clumped winds showed that 
wind clu mping may stro n gly af f ect the X-ray l i ne formation 
(see e.g. iFeldmeier et all 120031: lOskinova et all |2004 l2006t 
IOwocki & Cohen, 200^ 

The radiative transfer model s for n on-LTE stellar- 
winds, Uke CMFGEN (.HilUer & Milleil Il998h. PoWR 



(Haman n&Grafeneii 12004 and R\STWIND (iPuls et al.L 
|2005D, accaunt for wind inhomogeneities in an approximative 
way. An adjustable parameter, like the "clumping factor" D, 
defines the enhancement of the density inside clumps compared 
to a smooth model with the same mass-loss rate. Clumps are 
assumed to be optically thin at all frequencies, and the velocity 
field is monotonic. The empirical mass-loss rates derived under 
these assumptions (the so-called microclumping approach) 
differ from those obtained when different diagnostics are used. 
It was shown that the mass-loss rates derived from -based 
diagnostics (i.e. recombination lines such as Ha, IR, and radio 
emission lines) have to be reduced by a factor of Vd compared 
to the value s obtained under the a s sumpt i on of a smoo t h wind 
(see, e.g., iHamann & Koesterkel Il998t iBouret et all I2003L 
l2005h . On the other hand, diagnostics that depend linearly on 
the density (e.g., unsaturated UV resonan ce lines) sho u ld not 
be affected by optically thin clumping. Fulle rton et al.l (l2006l) 
found that mass-loss rates obtained from the Pv resonance 
fines are systematically smaller than those derived from Ha 
or radio free-free emission. Consequently, a reduction of the 
empirical mass-loss rates by up to a factor of 100 compared to 
the unclumped models was suggested. Another possibility to 
explain the discr epancies between p- and p^-b ased diagnostics 
was proposed bv IWaldron & Cassinelfil (l2010l) . They suggested 
that the XUV radiation near the He ii ionization edge originating 
in wind shocks may destroy the Pv ions and, consequently, 
Pv resonance line diagnostic may be explained without the 
need to d ecrease the mass-loss rates. However, it was shown by 
iKrticka & Kubat (2009) that X-rays fail to destroy the P v ions. 

In order to reconcile results obtained from different diagnos- 
tics, traditionally used assumptions hav e to be relaxed. The firs t 
attempt in this direction was made by lOskinova et al.l (l2007h . 
They studied resonance and recombination lines assuming sta- 
tistical properties of clumps and keeping a monotonic velocity 
field. They proposed that the discrepancies between mass-loss 
rates derived from recombination lines and from the P v reso- 
nance doublet can be solved by accounting for macroclumping. 
In this approach clumps can be of any optical thickness. They 
showed that accounting for macroclumping has a significant im- 
pact on the line formation process, manifested as reduction of 
the effective opacity of the medium leading to weaker lines for a 
given mass-loss rate. 

A further relaxation of the traditional assumptions was made 
byfOwocki (2008). He pointed out that for line transitions, the 
non-monotonic velocity field ("vorosity") can be important. He 
showed that a no n-monotonic veloci ty field may also reduce the 
effective opacitv. lZsargo et al] (|2008|) stressed that the void inter- 
clump medium (ICM) assumption also has to be relaxed, and 
that the ICM must be taken into account. Sundqvist et al. (2010) 
showed that the detailed density structure, the non-void ICM, 
and non-monotonic velocity field improve the line fits. They 



confirmed the findings of lOskinova et al.l (l2007h that the micro- 
clumping approximation is not adequate for UV resonance line 
for mation in t ypical OB-star winds. This is in agreement with 
IPrinia & Massal (l2010l) . who established spectroscopic evidence 
for optically thick clumps in the wind by measuring the ratios 
of the radial optical depths of the red and blue c omponents of 
the Si IV doublets of BO to B5 supergiants (see also Mass a et all 
2008). They showed that these ratios are spread between 1 and 
2. Since they differ from the predicted value of two (this value 
follows from the atomic constants for a smooth wind), this is a 
direct signature of optically thick clumping. 

In addition to the statistical approach, the Monte Carlo (MC) 
technique proved to be a very suitable method for studying 
clumps in hot star winds. lOskinova et al.l (l2004l2006h used this 
approach in 2D geometry to calculate X-ray line profiles. In a 
recent paper, Sundqvist et al. (2010) synthesized UV resonance 
lines from inhomogeneous pseudo-2D radiation-hydrodynamic 
wind models and 2-D stochastic wind models using MC radia- 



tive tr ansfer calculation. In their second paper, Sundqvi st et al 
(201 1) extended the wind model to pseudo-3-D. Muiires et aL 
(2011) also used a MC method combined with non-LTE model 
atmospheres to compute the effect of clumping and porosity on 
the momentum transfer from the radiation field to the wind. They 
parameterized clumping and porosity by heuristic prescriptions. 

Important steps has been made in the last five years in under- 
standing the line formation in clumped winds. However, all pre- 
vious works applied simplifications to the wind geometry. The 
structured stellar winds are essentially a 3-D problem, and a full 
description requires a 3-D radiative transfer. In this paper we 
treat for the first time the full problem with 3-D radiative trans- 
fer in the clumped wind. We accounting for: non-monotonic ve- 
locity, non-void ICM, and full 3-D without any limitation for 
geometry. We present the basic concept of the model and the 
method we developed. The main effects on the resonance lines 
(both singlets and doublets) including wind clumping in density 
and velocity as well as the effect of a non-void ICM are demon- 
strated. In a forthcoming publication we will apply our model to 
observations and derive mass-loss rates. 

In Sect. |2] we define the wind model (geometry, velocity, and 
opacity of the wind) and our parametrization of the clump prop- 
erties. The MC radiative transfer code is described in Sect. [3] 
Results of the model calculation are presented in Sect. H] In 
Sect. |5] we summarize our results and outline further work. 



2. The wind model 

We assumed a wind that may consist of a smooth and a clumped 
region. The clumped region comprises two density components: 
the tenuous ICM and the dense clumps. All distances in the wind 
are expressed in units of stellar radius. We introduce a parameter 
rci where the wind clumping sets on. The lower boundary of the 
wind r^in is set to the surface of the star {r^nin = !)■ The region 
I < r < Tci, represents the smooth wind region. The region < 
^ '"max, where rmax is the outer boundary of the wind, represents 
the clumped region (see Fig.[T]i. 

General method. In order to solve the radiative transfer 
through the clumped wind, we first generate a snapshot of the 
clumps' distribution. Using the MC approach (see Sect. [3]), we 
then follow the photons along their paths. The density and ve- 
locity of the wind can be arbitrarily defined in a 3-D space. The 
calculations are carried ou t in the comoving frame following the 
prescriptions bv lHamannI (fT98(3t) . 
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Fig. 1. A schematic view of the wind model. r^i„ is the lower 
boundary of the wind (the surface of the star), is the onset 
radius of the wind clumping, and r^ax is the outer boundary of 
the calculation. Vectors ri, rz, and k are described in Sect. 13.21 



We solve the radiative transfer in the dimensionless form. We 
introduce the dimensionless frequency in the observer's frame 
^obs^ measured relative to the line center in units of standard 
Doppler- width Avs, 



Wind velocity. The velocity field can be arbitrary. For simplic- 
ity, we assume that the velocity field is radial and for the smooth 
wind we adopt the standard yS-velocity law. 



(4) 



where u,- is the radial component of the velocity vector, Uco is the 
wind terminal velocity and b is chosen such that u^i„ < IQ^^Uoa- 
The derivative of the Eq. (|4|i gives radial velocity gradient 



, du, b 

V = = UrP , 

dr r(r — b) 

that is important for choosing an adequate integration step. 



(5) 



Wind opacity. Line opacity calculation in the moving medium 
requires that we take into account the Doppler shift, which 
makes the opacity anisotropic in the observers frame. We do not 
adopt the Sobolev approximation but solve the radiative transfer 
equation in the comoving frame. Since in this frame the fluid is 

at rest, the opacity is isotropic. 

Fo r the opacity we use the same parameterization as lHamaiiiil 
(fT980l) . 

Af(r)= 4^^(r)f,, (6) 



^obs 



(1) 



where is an arbitrary reference velocity and - /S.Vs{cIvq). 
During the calculations we strictly keep the local co-moving 
frame frequency jc*^™^ of the photon which means that the fre- 
quency of the followed photon at the certain coordinate point 
can be expressed using the scalar product of the vector of the lo- 
cal macroscopic velocity (d) and the unit vector in the direction 
of the photon's propagation (k). 



k ■ V. 



(2) 



Only when the photon exits the wind, x^"^^ must be transformed 
into the observer's frame by 



,,cmf 



k v. 



(3) 



Note that all velocities are dimensionless and measured in units 
of Vs. 

The code which we developed is gridless and does not re- 
quire any symmetry. Instead of using a predefined grid, we intro- 
duce an adaptive integration step Ar. After choosing Ar, the ve- 
locity and opacity are calculated along the photon's path. This al- 
lows us to account for arbitrary density and velocity inhomogen- 
ites (for a detailed description see Sect. 13.21 ). 

For the vectors we use a Cartesian coordinate system and in 
this coordinate system they are expressed by a set of coordinates 
(x, y, z). 

Basic assumptions. In order to study the basic effects of clump- 
ing on the resonance line formation (both singlets and doublets), 
we adopt a core-halo model. Only the line opacity is taken into 
account while the continuum opacity is neglected in the wind. 
We consider only pure scattering and assume complete redistri- 
bution. Only Doppler broadening is considered. 



where ;ifo is a free parameter of the model that corresponds to the 
line strength and it is proportional to the mass-loss rate and abun- 
dance of the absorbing ion, q(r) reflects the depth-dependent de- 
gree of the ionization of the absorbing ion and for simplicity it 
is chosen to be q{r) = 1 here, i.e. constant ionization condition. 
The absorption profile (px for the singlet lines is assumed to be 
Gaussian 

cf>.^-^e-''\ (7) 

where x is the dimensionless frequency given by Eq. ([TJ. We 
imply that the reference velocity is the Doppler-broadening 
velocity ud (i.e. - ud) that includes contribution from both 
thermal broadening and microturbulence and it is constant over 
the whole wind. 

In this parametric formalism, the smooth wind opacity xif) 
is proportional to the smooth wind density Psw('")- In our calcu- 
lations, we do not directly calculate wind density, but the depth- 
dependent line opacity ;if(r), which is enough for calculation of 
the emergent radiation from the wind. 



2. 1 . Description of clumping 
2.1.1. Clump properties 

Clumps are gaseous regions in the wind with higher density than 
their surroundings, and, possibly, also with velocities different 
from the smooth wind. 

We allow for arbitrary optical depth of clumps. The clumps 
can be optically thick in the cores of a resonance lines, while they 
may remain optically thin at all other frequencies. Clumps are 
assumed to be stochastically distributed within the stellar wind. 
The average clump separation L(r), measured between their cen- 
ters, is variable with the distance r from the star. For simplic- 
ity, we assume that each clump has the spherical shape with the 
volume {4-7r/3)P, where I is the radius of the particular clump. 
Clump radius varies with the distance from the star, / - l(r). The 
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Fig. 2. 2-D projection of an example of a realisation of our stochastic 3-D wind model. Left: Clumps (blue circles) with the different 
size distributed around the star (red filled circle). Right: Distribution of the clumps; the filled red sphere in the center represent the 
star, while the blue dots represent the positions of the clump centers. 



density pci inside the clump at a distance r from the center is as- 
sumed to be by a factor D higher than the smooth wind density 
Psw at the same radius. 



(8) 



and D > 1 . The factor D (density contrast) is the first free pa- 
rameter of our model and for simplicity it is assumed to be depth 
independent. Using L and I, the factor D can be expressed as 



D = 



3 



l\r) 



(9) 



In our model clumps are assumed to be preserved entities. We 
assume that clumps are accelerated radially and neither split nor 
merge. The consequence of this is that the number of clumps per 
unit volume «ci has to satisfy the equation of the continuity. 



1 



(10) 



Then for the average clump separation, holds L 
Eq. ( [Tol l, this can be expressed as 



n^'^. With 



(11) 



where Lq is a second free parameter of our model that represents 
THE typical clump separation in units of the stellar radius, and 
\D{r) - Vrlvac is the velocity in the units of the terminal speed. 
Similar as in the Eq. (fTTI) . the clump radius can be written as 



/(r) = /o^/r2 w{r), 



where 



AnD 



(12) 



(13) 



which follows from Eqs. ^ and (fTTT l. With this procedure we 
create clumps whose size varies with the radial distance (see 
Fig. 13. 



The volume filling factor fy is defined as the ratio of the 
total volume that clumps occupy (Vd) to the whole volume y„ 
in which the clumps are distributed (between and rmax), 



„ Vcl 1=1 

JV - — - 



^ y max ' q\J 



(14) 



where Vd, = {Ajil3>)l] is the volume of /-th clump and A^ci is the 

total number of the clumps. 

For the case of void ICM, the entire mass of the wind m„ = 
Psv/ir) dV - {psw) Vv/ in the clumps can be written as (using 

Eq.® 

Nd Nd Nd Nd 

»Jw = ^ '«cl; = ^ Pel/ Vcl; = D^p'^^VcU = D <psw> ^ Vcl; 



(=1 



(=1 



/=1 



/=1 



(15) 

where m^u and pd, are mass and density of the /-th clump, re- 
spectively, and p'^^ is the value of the smooth wind density at the 
location of the /-th clump. If the ICM is void. 



(16) 



2.1.2. Inter-clump medium 



The ICM density pic is assumed to be reduced with respect to 
the smooth wind density psw (the density of the wind without 
any clumping) by a factor d (third free parameter of our model, 
assumed to be radius independent). 



Pic(r) = psw(r); 0<<i<l. 



(17) 



For the case of non-void ICM, the total mass of the wind is 
distributed between clumps and ICM. In this case, the total mass 
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Fig. 3. Velocity structure of the wind model along the considered 
photon path for the case of non-monotonic velocity distribution 
inside clumps. 



where is the absolute position of the center of the i-th clump, 
v/iir^) is the velocity determined according to Eq. ^ at the 
position r^. The velocity dispersion i'dis('") - ntvpif) where m 
(0 < m < 1) is the velocity deviation parameter (fourth free pa- 
rameter of our model). Eq. ( l20b introduces a negative velocity 
gradient inside clumps, while the center of clump has the ve- 
locity as the smooth wind. The velocity structure of our model 
along one particular photon path is shown in Fig. [3] 



3. Monte Carlo radiative transfer 

The MC approach is the most promising way for treatment of 
radiation transfer in a clumped stellar wind. "Classical" solu- 
tion methods for the radiative transfer equation, like the Feautrier 
scheme or short characteristics methods, become extremely time 
and memory consuming and need very sophisticated solution 
schemes when we go to more spatial dimensions than one (e.g., 
lKorcakova& KiiMll200l 12001 iLobel & BlommeLl200§ ). This 
is not the case of MC methods, where the extension to 2-D or 3- 
D is relatively simple. 



of the wind - (psw) can be expressed as 

Nd A'd 

mw = mic-l-^ m^i = d{p^^){V^-Vd) + D {p^^)^ Vd/. (18) 
Then from the Eqs. ( fT4l i and ( fTST i follows 

2.1.3. Clump distribution 

The distance r, of the /-th clump from the stellar center is cho- 
sen randomly. By e mploying the von Neumann rejection method 
dPress et al.L 1 1992b with the inverse of the velocity law l/y^ as 
the probability density distribution function. This reflects the 
equation of the continuity for the number density of clumps. 
Consequently, more clumps are concentrated close to the star 
(see Fig.|2]i. 

Clumps are distributed uniformly in cos 0, and (^,, and 6, and 
<fi are randomly chosen with cos 9i - /v; - and if, - 2n^a 
i^n and ^,2 are two diff'erent random numbers). For given r,, the 
radius of each clump /, is determined according to Eq. (fTST i. We 
do not allow the clumps to overlap. 

2.1.4. Inliomogeneous velocity in clumps 

In our wind model the velocity of the smooth and the 
inter-clump medium is assumed to be monotonic, Eq. (HJi. 
However, the velocity inside the clumps is allowed to devi- 
ate from the monotonic wind. Consequently, a negative ve- 
locity gradient can appear (see more in the Sect. O. This as- 
sumption is base d on the prediction of hydrodynami c wind 
simulations (e.g. ,O wocki et al.l 119881; iFeldmeier et"^ Il997bt 
iRunacres & Owock i, 2002), which show that over-dense regions 
inside the wind are slower and that ICM has approximately the 
same velocity as the smooth wind. 

The velocity inside the /-th clump can be approximated as 

v(r) = vi^ir'r) - va,,(r) (20) 



3.1. Creation of a photon 

A photon is released from the lower boundary of the wind (sur- 
face of the star) and then follow through the wind until they reach 
the outer boundary of the wind or scatters back into the photo- 
sphere. The photons are uniformly distributed in cos 6 and tp over 
the whole surface area of the star. The photons from the lower 
boundary are released only upwards uniformly in (p and with a 
distribution function oc yuc/yu in fi. The angular distribution func- 
tion for photons emitted by the photosphere follows from the 
definition of the flux (see, e.g., Lucv, 1983). The initial unit vec- 
tor k of the photon propagation from the surface of the star is 
randomly chosen with cos^i - fi^ = VST and tpk = 27r^t2 i^ki 
and are two different random numbers). 

In our MC calculations, frequencies of newly created pho- 
tons are det ermined from t he interval defined using the ratio 
Vco/vd (as in lHamarinIll98lh . Initial frequencies of photons have 
values from the interval (jCmin, JCmax) where 

— _ _ ■ — 4. ('91^ 

-'^min — .^band i -'^max — ~i" .^band ■ \^ ^) 

Vd Vo 

where Xband is the line width and it is usually assumed t o be about 
4.5 Doppler units (similarly as in Mih alas et al. 

We divide the whole frequency interval into A^bin subinter- 
vals of equal length, and the same number of the photons are 
released in each subinterval. Each photon obtains the frequency 
xf'^ir^ia), randomly chosen from the frequency interval of n-th 
(« = 1 , . . . , A^bin) bin, given in the observer's frame as 

xf'^r^m) = (n - ^)xbin + Xmin, (22) 

where Xbin = (jCmax - JCmin) /A^bin IS the width of the frequency 
bin (the same for all bins) and ^ is a random number. Then the 
frequency is transformed to the comoving-frame frequency x™^. 

3.2. Optical deptii calculation 

After its creation, each photon obtains an information how far 
it is allowed to travel before it undergoes interaction with some 
particle. This distance is d etermined by a randorn l y chosen op - 
tical depth = -In^ (see lAverv & Housel [T968t fcaroff et all 
11972) . Then the actual optical depth, which photon passes on its 
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Fig. 4. The path of one particular photon inside a reahzation of 
our clumped wind with the adaptive integration step. The bigger 
green sphere represents the outer boundary while the red sphere 
in the center represents the lower boundary of the wind. Smaller 
blue balls represent the clumps and black dots denote integration 
steps. 



travel, is calculated by summing opacity contribution along its 
path, and it is checked with r^. 

In order not to miss any inhomogeneities inside the wind in 
either velocity or density, we use adaptive integration steps. The 
photon path is defined by the unit vector of photon propagation 
k and the photon is followed along its path by stepwise adding a 
variable integration step Ar. If the current photon position is ri , 
then the vector of the next photon position r2 is r2 = ri + Ar k. 
After each integration step, the coordinates of the current photon 
position are recalculated. When the photon is followed through 
the smooth part of the wind, Ar is set not to be smaller than 
0.1 of the radial length of the smooth wind region and not to 
be larger than Ur/u', i.e., 0.1 (r^ - rmin) < Ar < Ur/u', u' is the 
velocity derivative - see Eq. (|5]l. When the photon is followed 
inside the ICM, then Ar = 0.3 /(ra), where /(ra) is the radius of 
the clump that would have at r2 (Eq. fT2l l. For the photon path 
followed inside the clump, Ar is not smaller than 0. 1 l(r2) and 
not larger than 0.3 /u', i.e., 0.1 /(r2) < Ar < 0.3 /v'. This adaptive 
integration step technique ensures that opacity calculation along 
the photon path is done properly without a chance to skip any 
inhomogeneity inside the wind and that it is fast enough. The 
illustration of the adaptive integration step is shown in Fig.|4] 

For the optical depth calculation, it is necessary to know the 
opacity. According to Eq. for the opacity calculation at the 
current position we need to know i>r and the Doppler shifted fre- 
quency of the photon at every position of the integration process. 
When the photon is traced inside the smooth or inter-clump re- 
gions, the velocity is calculated according to Eq. (|4|. When the 
photon is traced inside the i-th clump and the option of the inho- 
mogeneous velocity is turned on (Sect. 12.1.4b . then the velocity 
is calculated according to Eq. ( l20b . 

Along the photon path the Doppler shifted comoving-frame 
frequency of the photon is calculated. 



The optical depth between ri and for the actual integration 
step j is 




and it is calculated using the trapezoidal rule. Along the whole 
path, the optical depth is accumulated (r = 2 j tj, where the J 
is the total number of the integration steps made before scatter- 
ing happens), and when the total optical depth t > r^, then the 
condition for the line scattering is fulfilled. 

After the scattering, the photon obtains a new direction (9, (p), 
chosen randomly for the case of isotropic scattering as cos = 
2^-1 and (p = In^ and a new optical depth = - In ^ is ran- 
domly chosen again. For determination of 9, 0, and three dif- 
ferent random numbers are used. Assuming complete redistribu- 
tion the photon frequency after scattering is calculated randomly 
from Gaussian distribution with the mean zero and the standard 
deviation unity (Box & Muller, 1958). 

In order to reproduce the emergent flux from the wind, all 
photons that reach the edge of the wind are collected and put 
into a proper frequency bin. Each of the profiles shown in this pa- 
per has been calculated for one random configuration of clumps. 
This would correspond to an observation with a short exposure 
time, compared to the dynamical time scale. But on the other 
hand, we count the emergent photons irrespective of their direc- 
tion, while the distant observer sees the 3-D clump configuration 
from one specific direction. By this averaging over all directions, 
we effectively obtain a mean emergent profile, similar to what 
we would get from averaging over many different random clump 
configurations but seen from one specific direction. 

4. Results of the 3-D wind model calculation 

In this section we show basic effects on the line profiles for both 
single resonance lines and resonance doublets varying different 
model parameters. Depending on which effect we want to study, 
some parameters are kept fixed, while others are varied within 
considered ranges (see Table [TJ. The basic model parameters 
(left part of Table [U are chosen to represent a typical O type 
star The varied model parameters (the right part of Table [l]) 
have their default value except the selected one whose effect we 
aim to study. The effects on the line profiles are shown for weak 
(Xo = 1), intermediate (yo = 10), and strong (xq - 100) lines. If 
the variation of some model parameter shows a similar effect on 
all three types of lines, we show that effect only for the strong 
line case. 

All calculations are performed using 10^ photons distributed 
over 100 frequency bins, resulting in a S/N ratio of about 30 per 
bin due to the Poisson statistics. 

4. 1 . The effects of the macroclumping 

The total number of clumps in one snapshot of clump distribu- 
tion can be derived using the parameter Lq from 

'max 

A^ci = j n^xir)AnP-dr. (24) 

For instance, if 0.2 < Lq ^ 0.5, then the total number of clumps 
within the clumped region is 10^ > A^ci <: 10^. A higher number 
of clumps (i.e smaller values of Lq) produces a less porous wind. 
If Lo — > the whole wind is filled with many little clumps and it 
resembles a smooth wind. 
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Fig. 5. The effects of the macroclumping (left) and the non-void ICM (right) on the weak (xo = 1, the upper panels), intermediate 
(Xq = 10, the middle panels), and strong (yo = 100, the lower panels) lines. Left: The black dashed lines represent a smooth wind 
(Lo 0) and other lines are calculated for a different clump separation parameter Lq as given in the panels. Other model parameters 
have their default value (Table[T]). Right: The black dashed lines represent a smooth wind and other lines are calculated for different 
values of the ICM density parameter d as given in the panels. Other model parameters have their default value (Table[T]). 



In this subsection, all model parameters are set to their de- 
fault values except of xa and Lo, which are varied as given in 
Table [T] The main macroclumping effect on the line profile is 
the reduction of the line strength compared to the smooth wind. 
For this case, where we assumed that clumping starts from the 
surface of the star (rd - 0), the variation of the Lq parameter has 



a strong influence on all three types of lines (see the left part of 
Fig. |5]l. This is because the clumps are optically thick for these 
lines. However, if we assume that clumping starts at r^\ - 1.3, 
the clumps are optically thin for the weak lines there while they 
remind optically thick for the intermediate and strong line case 
(see lSurlanetaLl 120121) . 
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Table 1. The model parameters to study the influence of their variation on the resonance line profiles. 
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Fig. 6. The effect of variation of the clumping factor D for the 
case of a strong line (xo - 100). The black dashed line rep- 
resents a smooth wind, while the other lines are calculated for 
different D as given in the figure. Other model parameters have 
their default value (Table [T]i. 
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Fig. 7. The effect of variation of the onset of the clumping for 
the case of strong line (xo - 100). The dashed black line repre- 
sents the smooth wind (r^ - ^max) and other lines are calculated 
for different r^i as given in the figure. Other model parameters 
have their default value (Table [T]l. 



According to Eqs. (fTSl i and ( fT3] ), the radius of the clumps 
depends on Lq. For very small Lq a huge number of clumps with 
small radii are created, and the clump contribution to the line 
opacity is almost the same as in case of a smooth wind. In this 
case there are not too many "holes" between the clumps, and 
the photons cannot escape from the wind easily. But for higher 
values of Lq, less clumps with larger radii exist. Consequently, 
there are more "holes" in the wind, through which photons may 
freely propagate. This leads to lower absorption and weaker line. 

Clumping lowers the effective opacity. This effect is weaker 
for the outer parts of the wind, because the individual clumps be- 
come optically thin there. This causes the bump in the blue part 
of the line. On the other hand, a significant part of the wind has 
approximately the same wind velocity close to the wind termi- 
nal velocity. Consequently, despite lower effective opacity there 
is still enough matter to absorb. This causes the absorption dip 
near Voa especially for lower Lq. 

A similar effect of attenuation of the line strength can be 
obtained by varying the clumping factor D. In Fig. |6] we show 
the effect only for the strong line (xq - 100). The model pa- 
rameters are again set to their default values except of D that 
is varied. When enhancing the density inside clumps compared 
to the smooth wind density (i.e. increasing the parameter D), the 
porosity effect is more pronounced and the line becomes weaker. 



4.2. Effects of the non-void inter clump medium 

In this subsection we study the changes caused by variations 
of the ICM density parameter d and xo, while the other model 
parameters are set to their default values. The space between 
clumps is filled with some amount of matter with density de- 
fined by d. This matter fills the density "holes" between the 
clumps, and photons can also be scattered there. This manifests 
as strengthening of both the absorption and emission parts of 
the line profiles with respect to the model with void ICM (see 
the right part of Fig. |5]l. This effect is most pronounced for the 
strong lines, where even for a small d the ICM contributes a 
lot to the absorption coefficient and makes the lines more satu- 
rated and the emission parts stronger (see the right lower panel 
in Fig.jSj. For the same value of the parameter d, the strong lines 
are saturated more than the intermediate lines. A proper choice 
of d can saturate the strong lines but still keep the intermediate 
lines unsaturated (e.g. see lines for = 0.1 in the right middle 
and lower panels of Fig.|5]l. 

Assuming that the ICM is not void has an influence on the 
strength of the line as a whole, but it has a particular effect on the 
center of the line, where a small absorption dip appears. Because 
more clumps are concentrated close to the star, the inner part of 
the wind together with the ICM contribute to the line opacity 
more than the outer part of the wind. 
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4.3. Effects of the onset of clumping 

In order to show the effect of the radius where clumping starts, 
we now fix all parameters of the model at their default values 
except the onset radius of clumping r^. The main effect of the r^i 
variation is the appearance of a strong absorption near the line 
center (see Fig. |7]), which is due to the smooth part of the wind 
inside r^i at low expansion velocity. If clumping starts higher 
in the wind, the absorption near the line center is broader and 
also partly attenuates the emission from the back hemisphere of 
the wind. Therefore, this central absorption dip can be used as 
diagnostic for the onset of clumping. 

It is interesting to note that even a very little smooth part of 
the wind (between r = 1 and r - 1.01) causes significant central 
absorption, an effect which should be observable. As it can be 
seen from the Fig.]?] only when rd - 1 the absorption near the 
line center disappears. When the clumped part of the wind is 
larger (smaller value of the r^i), the reduction of the line strength 
is more pronounced. By setting r^i = r^ax, the smooth wind is 
reproduced. 

4.4. Effects of velocity dispersion inside clumps 

The absorption of an individual clump can be broadened 
by stochastic (thermal or micro-turbulent) motions inside the 
clump, but also by an additional velocity gradient in the clump 
as predicted by hydrodynamic simulations. These kinds of 
line broadening are described in our model by the Doppler- 
broadening velocity (uu), and by the velocity dispersion inside 
the clumps (wdis), respectively (ses Sect. 12. 1.4b . 

In Fig. [8] we show the effect of the vu variation. The nor- 
malization of the frequency-integrated absorption coefficient is 
maintained by a compensating change ofxo- For the higher value 
of V£, the line profile is broader, and absorption and emission are 
stronger. When decreases, the macroclumping effect becomes 
more pronounced because the clumps are optically thicker in the 
line center, which leads to a smaller effective optical depth due 
to a larger macroclumping effect. 

To study the effect of the velocity dispersion inside clumps, 
we fix uu = 20kms"' and vary m and ;t'o (Fig- HI- If the ve- 
locity dispersion inside the clumps is higher (i.e. when increas- 
ing the parameter m), the gaps in the velocity field are smaller 
(more velocities overlap) and the probability of photon escape 
is lower. This leads to some absorption at velocities higher than 
Voa. However, in the shown example the velocity gradient inside 
the clumps does not differ much from the velocity gradient of the 
smooth wind except of its sign, thus leaving the optical depth of 
the individual clumps is roughly the same. Therefore, the "voros- 
ity" has only little effect on the total line strength. 

When the velocity dispersion is accounted for, the absorption 
extends to velocities higher than the terminal velocity Uco- This 
eff ect is important for the derivation of from observations (i.e. 
see iPrinia et all Il990h . and may provide an effective diagnostic 
for the existence of the velocity dispersion inside the clumps. 

4.5. Doublet calculation 

For the calculation of doublets, the profile <p^ in Eq. (|6]l is as- 
sumed to have the form 

0, = — (e-(-'^+''»p/2)' + p e-(-»-*=p/2)')^ (25) 

where lisep is the separation between doublet components. The 
ratio of line opacities has a fixed value p - Xi/X2, which fol- 
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Fig. 8. The effect of different Doppler broadening on the line 
profiles represented by three values of and XQ given in 
the figure. Other model parameters have their default value 
(Tableffli. 

lows from the atomic line strengths. The zero point frequency is 
located in the middle between the components. 

The frequency of the photon after scattering is chosen ran- 
domly with a Gaussian distribution in the same way as for sin- 
glets. However, if the co-moving frame frequency of the photon 
before scattering indicates that it belongs to the blue component ( 
j^cmf J. Q-)^ assume that the photon is redistributed only within 
the blue component of the doublet, and vice versa. 

In Fig. [To] we show one example of a strong doublet with 
(isep - 2000 km s"' and;i'o = 500. We demonstrate different ef- 
fects on the line profile by changing the properties of the clumps. 
These effects are analogues to the single-line case. The main ef- 
fect of macroclumping is the reduction of the line strength. If the 
velocity dispersion inside clumps is taken into account, there is 
some absorption at velocities higher than Uoo. Only for the case 
of non-void ICM it is possible to saturate the line. 

5. Summary 

In this paper we present a full 3-D inhomogenous stellar wind 
model and solve the radiative transfer to model the resonance 
lines. To our knowledge this is the first work where the prob- 
lem of resonance line formation in stellar winds is solved in full 
3-D, while previous work was restricted to 2-D or pseudo-3-D 
geometries only. 

The radiative transfer is also calculated for the formation of 
doublets. This is important, because UV resonance doublets are 
a key diagnostic for stellar winds. Our work demonstrates the 
role which stellar wind clumping plays in the formation of spec- 
tral lines, and how it affects the empirical wind diagnostics. 

The method we develop here is very flexible, as it is capable 
to account for different 3-D shapes of clumps, for arbitrary 3- 
D velocity fields, and for ICM the same time. Below we briefly 
summarize the main results of our models, deferring the appli- 
cation to observed spectra to forthcoming work (Surlan et al., in 
prep.). 

• Allowing for clumps of any optical depths (macroclump- 
ing) causes a reduction of the effective opacity in the lines. 
Consequently, there is less absorption in the wind. Since this 
reduction is weaker for the outer parts of the wind, the line 
profiles show an absorption dip near Doo. 
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Fig. 9. Effects of the velocity dispersion inside clumps on the 
line profile. Upper panel: weak lines (xq - 2.5), middle 
panel: intermediate lines (xq - 25), lower panel: strong lines 
iXo = 250). The black dashed lines represent the smooth wind, 
the green dashed lines (poro) represent pure porous wind, the 
other lines (poro+voro) represent the porous wind with non- 
monotonic velocity described with m as given in the panels. 
vu = 20kms ' and other model parameters have their default 
value (Table [TJ. 
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Fig. 10. Effects of the macroclumping on the strong doublet line 
profiles (xq - 500) including non-void ICM and velocity disper- 
sion inside clumps. The black dashed line represents the smooth 
wind, the green dashed line (poro) represents the macroclump- 
ing effect, the full red Une (poro-Hvoro) represents the macro- 
clumping effect including non-monotonic velocity described in- 
side clumps with m - 0.2, and the orange dash dotted line 
(poro-fvoro+icm) represents the macroclumping effect includ- 
ing non-monotonic velocity with m - 0.2, and non-void ICM 
with d - 0.05. Other model parameters have their default value 
(Table©. 



For a given clumping factor D, the key model parameter af- 
fecting the effective opacity is Lq, the clump separation pa- 
rameter. The opacity reduction is largest for the largest Lo 
(Fig. |5]l. Therefore, the mass-loss rate empirically obtained 
by fitting the UV resonance lines becomes larger when the 
macroclumping effect is taken into account. 
The onset of clumping, r^\, affects the line shape: the closer 
to the stellar surface clumping starts, the more pronounced is 
the absorption dip at the line center (Fig.|7]l. This absorption 



dip may provide an effective a diagnostic for the onset 
of clumping. 

• The line saturation is strongly affected by the ICM. A non- 
void ICM is required to reproduce the saturated lines simul- 
taneously with non-saturated lines (Fig.|5j. 

• When accounting for a velocity dispersion within the 
clumps, added to the mean velocity law, the absorption ex- 
tends to a larger blue-shift than corresponding to Vco- This 
effect has to be taken into account when deriving v^o from 
observations. 

• In any clumped wind, non-monotonic velocities will always 
appear together with the density inhomogeneities. Therefore, 
their combined effect must be taken into account for the Une 
formation modeling. 

• In case of resonance doublets, the clumping effects are ana- 
logues to the case of single lines. 

The main conclusion of our work is that in a realistic 3-D 
wind with density inhomogeneities and non-stationary veloc- 
ity, the P-Cygni profiles from resonance lines are different from 
those from smooth and stationary 3-D winds. Any mass-loss di- 
agnostics which do not account for wind clumping must under- 
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estimate the actual mass-loss rates. This can explain the reported 
discrepancies between the mass-loss rates obtained from p- and 
p^-based diagnostics, respectively. Using our general description 
of clumping presented in this work, it will be possible to deter- 
mine improved values of the stellar mass-loss rates. 
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